Correlated band structure of electron-doped cuprate materials 
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i We present a numerical study of the doping dependence of the spectral function of the n-type 

. cuprates. Using a variational cluster-perturbation theory approach based upon the self-energy- 

■ functional theory, the spectral function of the electron-doped two-dimensional Hubbard model is 

£Nj ' calculated. The model includes the next-nearest neighbor electronic hopping amplitude t' and a 

fixed on-site interaction U = 8t at half filling and doping levels ranging from x = 0.077 to x — 0.20. 
Our results support the fact that a comprehensive description of the single-particle spectrum of 
electron-doped cuprates requires a proper treatment of strong electronic correlations. In contrast to 
previous weak-coupling approaches, we obtain a consistent description of the ARPES experiments 
without the need to introduce a doping-dependent on-site interaction U. 
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INTRODUCTION 



Angular resolved photo-emission spectroscopy (ARPES) has greatly contributed to our current understanding of 
materials with strong electron correlations, in particular high-temperature superconductors (HTSC). Up to a few 
' years ago, ARPES investigations have been concentrated on hole-doped HTSC materials. Since ARPES probes 
(—] . the part of the spectral function that is occupied by electrons, only the region below the insulating gap can be 
investigated in hole-doped cuprates. Although the observation of the unoccupied parts of the spectral function is in 
principle possible via inverse photoemission, the process is highly involved and does not yield the same resolution 
as direct photoemission. An opportunity for a more comprehensive study of the doping dependence of the spectral 
— i . function is offered by electron-doped cuprates. In these materials, not only the excitations below the Fermi level 
in the lower Hubbard band, but also those below the Fermi level in the upper Hubbard band can be studied and 
thus a large part of the important low-energy excitations is covered. For this reason, the investigation of such n-type 
cuprates by ARPES provides a large amount of new information on the quasiparticle dynamics of HTSC cuprates. 
■^j- . Recently, an ARPES study of the dopingdependence of the electron-doped cuprate Nd2_ a; Ce a ;CuOCl4±5 (NCCO) 
has been carried out by the Stanford group |l|,[2j. In these measurements, the low-energy excitations of Nd2CuOCl4±a 
(NCO) were shown to essentially coincide with the ones of typical undoped parent compounds of hole-doped materials 
such as Sr2Cu02Cl2 and Ca2Cu02Cl2, thus demonstrating the universality of the electronic structure of the (single 
layer) undoped cuprates. Upon electron doping, a remarkable Fermi surface (FS) evolution was found: In the heavily 
underdoped region the low-energy spectral weight is limited to an area close around k = (7T,0). This has been 
interpreted as the formation of electron pockets. With increasing doping level, these pockets are connected by FS 
' patches and finally form a large LDA-like FS closed around k = (ir, tt). 

The spectral function and FS data presented in pj gave rise to several theoretical interpretations, which also include 
the idea of a collapse of the Mott gap due to strong suppression of the local Coulomb repulsion upon doping^]. This 
conclusion was based upon mean-field calculations which employ a fitting of the on-site repulsion U e g of the Hubbard 
model as a function of doping, to the value of the antiferromagnetic order parameter obtained from experiments at 
^ , each doping level. 

In the present paper, using both standard cluster-perturbation theory (CPT) Q| 0. 0] and, in particular, also 
a variationally improved version (V-CPT) 0> HI 0> E3 > which is based on the self-energy functional theory we 
calculate the spectral function of the two-dimensional one-band Hubbard model. This model is taken with next-nearest 
neighbor hopping amplitude t' — — 0.35i and fixed on-site interaction U = 8t at half filling and doping levels ranging 
from x — 0.077 to x = 0.20. Our numerical results show that the salient features of the recent ARPES experiments 
for electron-doped cuprates can be reproduced with one-and-the-same Hubbard model without the necessity to resort 
to any change of the [/-values as a function of doping as used in previous theoretical studies. Our challenge here is to 
reproduce the global (i.e. n- and p-doped) phase diagram by one universal choice of the model parameters, starting 
from a picture of a doped Mott-Hubbard insulator. 



2 



MODEL 

As a generic model for the HTSC compounds we use the one-band Hubbard model 01 H\b, i-c. 

Hu, = -t (4a c j<r + h - c ) +U^2n i]nil . (1) 

(id) » 

Here, c| CT (c^) creates (annihilates) an electron on site i with spin a, (...) denotes nearest neighbors and U is the on-site 
part of the Coulomb repulsion. Although the t— U Hubbard model at low temperature develops a quasiparticle band of 
the appropriate width |l3|,[l4|, the dispersion shows a near degeneracy between the k = (tt, 0) and the k = (tt/2, 7t/2) 
points. From ARPES experiments, however, we know that the quasiparticle peak at k = (tt, 0) is shifted to higher 
binding energies. Actually, the dispersion of the quasiparticle peak shows two parabola with lowest binding energy at 
k = (tt/2,0) and k = (it, tt/2). It is, thus, indispensable to add at least one additional term ^3,0], which is taken 
here to be the hopping between next-nearest neighbors (((...))), i.e. 

- *' E ( C Ua + h.c) . (2) 

((id)) 

Even longer-range hopping (t") elements have been proposed to achieve consistency with experiment 0|. However, 
for the purpose of a qualitative analysis, it is sufficient to lift the degeneracy between k = (it, 0) and k = (tt/2, tt/2) 
and, thus, create the indirect gap as observed in experiments. 



NUMERICAL TECHNIQUE 



Despite the considerable simplification arising from the use of an effective single-band model, the calculation of 
the spectral function of the Hubbard model still remains a difficult task. The calculation of this quantity by exact 
diagonalisation is only possible for a small lattice of up to 4 x 4 sites, provided periodic boundary conditions are 
used. Larger lattice sizes can only be calculated by stochastic methods, such as the quantum Monte Carlo (QMC) 
technique |l8| or the density matrix renormalization group algorithm (DMRG) 0] . While these techniques certainly 
represent powerful approaches to strongly-correlated electron systems, they are known to show disadvantages when 
applied to the present problem. In case of QMC, doping and low temperatures lead to the well-known sign problem, 
i.e. the computation time increases exponentially with T and system size ^l|2(|. DMRG, in contrast, is a ground 
state technique successfully applied for 1-D and ladder systems, but displays convergence problems when applied to 
two-dimensional systems. 

Recently, a strong-coupling perturbation theory has been developed for which the infinite lattice is subdivided into 
sufficiently small clusters that can be treated exactly, followed by an infinite-lattice expansion in powers of the hopping 
between the clusters The expansion in the intercluster hopping can be formally carried out up to arbitrary 

order following the diagrammatic method of Refs. 

mini!. 

The lowest order of this strong-coupling expansion in 
the inter-cluster hopping has been termed "cluster perturbation theory" (CPT). 
The CPT Green's function is given by 

Goo = Go + GqTGo + GqTGqTGo . . . = [G^ 1 - T] _1 , (3) 

where G^ is the Green's function of the infinite system, Go the cluster Green's function and T the inter-cluster 
hopping. All quantities are matrices with indices referring to the particular cluster and to the sites within that 
cluster. CPT can be considered as a systematic approach with respect to the cluster size, i.e. it becomes exact in the 
limit N c — > oo, where N c is the number of sites within a cluster. In addition, CPT provides approximate results for 
an infinitely large system: One of the advantages of this method is that the CPT Green's function is defined for any 
wave vector k in the Brillouin zone, contrary to common "direct" cluster methods, like QMC or ED, for which only 
a few momenta are available. 

CPT results for static quantities as well as for the single-particle spectral function have been shown to agree very 
well with different exact analytical and numerical results |5|,|6|- On the other hand, there is also a serious disadvantage 
of CPT at this level: Namely, the method does not contain any self-consistent procedure which implies that symmetry- 
broken phases cannot be studied. We have recently proposed a variational approach (V-CPT) to this problem, which 
is based on the self-energy-functional approach (SFA). This method is explained in detail elsewhere |7l l8l 1 1 lj . We use 
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FIG. 1: Spectral function of the t - t' - U Hubbard model with t' = -0.35t and U = 8t. Panel A: half filling. Panel B: 
overdoped system with x = 0.2. Panel C: Detailed dispersion of the quasiparticle band of panel A for t = 0.5eV. Symbols 
represent the experimentally determined dispersion for Sr2Cu02Cl2 (squares circles |2ql . diamonds [2^1 ). Both spectra 
where obtained by a CPT calculation using a 10 site cluster. 



V-CPT to calculate the Green's function of the half-filled Hubbard model with long-range antiferromagnetic order 
while plain CPT is used for the doped system. The extension to a doped system with d-wave superconductivity was 
first carried out in Ref. see also Ref. |Toj | . 



RESULTS AND DISCUSSION 



We consider the above single-band Hubbard model with nearest- t and next-nearest-neighbor t' hopping at zero 
temperature. Useful parameterizations of the t — t' — U Hubbard model can be taken from the literature [lfil E3 | . 
We choose t' = — 0.35i and U — 8t here, which yields a sufficiently accurate ratio for the Mott gap A rj At and for 
the width of the quasiparticle band W ~ It, and fits the experimental dispersion of the quasiparticle band. For a 
reasonable approximation to the full many-body problem, CPT and V-CPT calculations for relatively large clusters 
are required. In this work, we have calculated the spectral function for half filling and for doped systems with x = 0.07 
to 0.2. To achieve the smallest doping level, clusters consisting of 13 sites have been used. 

Figure n displays the spectral function A(k, ui) of the half-filled (x = 0.0, panel A) and overdoped t — t' — U Hubbard 
model (x = 0.2, panel B). The plots show A(k, cu) along the momenta r = (0, 0) — > X = (n, 0) — » M — (tt, tt) — > T = 
(0, 0) through the Brillouin zone. The half-filled system in panel A gives rise to a narrow quasiparticle band, roughly 
between u> — —St and lj = —2t. A more detailed plot is given in panel C. One notices the characteristic parabolic 
dispersion close to k = (7r/2,0V, k = (tt,ttJ2) and k = (ir/2,ir/2). Assuming t w 0.5eV, this dispersion is practically 
identical to the ARPES data [3, EE E3- . The indirect single particle gap between k = (tt, 0) and k = (tt/2,tt/2) 
is about At, which is the maximum value still compatible with experiments. 

At about uj = —3t a feature with maximal spectral weight close to k = (0, 0) can be observed. The same spectral 
feature was already observed in early QMC simulations of the single-band t — U Hubbard model [2^ , exact diagonal- 
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izations of the t-J model[30j and approximate perturbative methods such as the self-consistent Born approximation 
(SCBA). Its spectral weight was mostly perceived as coherent, i. e. corresponding to the "coherent" motion of a 
"spin-bag" like quasiparticle. This is also supported by more recent QMC simulations^, V-CPT calculations of 
the t — U Hubbard model, and analytical considerations j3l| . 

We now discuss the spectral function for the overdoped (x = 0.2) system, which is plotted in panel B. Here, one 
finds a metallic quasiparticle band with a flat dispersion just below the Fermi level at k = (n, 0). The band crosses 
the Fermi level close to k = (n/2, 7r/2) and k = (tt,tt/2) and therefore creates a large Fermi surface closed around 
k = (-7T, it) in the Brillouin zone. The quasiparticle band shows almost the same dispersion as the tight-binding 
(U = 0) model with the same parameterization. Deep below the Fermi level, between uj = —Qt and oj = — lOt, we can 
see the traces of the lower Hubbard band observed at half filling, i.e. the area between uj = — 4/j and uj ~ —6t in panel 
A. 
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FIG. 2: Details of the spectral function of the t — t' — U Hubbard model with t' = — 0.35t and U — 8t restricted to a narrower 
region around the Fermi level. Panel A: x=0.077, obtained from a 13 site cluster. Panel B: x=0.091, obtained from a 11 site 
cluster. Panel C: x=0.200, obtained from a 10 site cluster. 
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FIG. 3: Spectral weight at the Fermi energy, whose maximum is an indication of the Fermi-surface contour for different doping 
levels. The plots A through C are obtained from ARPES experiments (taken from |2|]), while D through F correspond to our 
CPT results, whose spectral functions are plotted in Fig. [5] 

In Fig. |21 we plot a detail of the spectral function of the Hubbard model in a more restricted region around the 
Fermi level for three doping values, namely x = 0.077 (panel A), x = 0.091 (panel B), and x = 0.200 (panel C). For 
x = 0.077, the Fermi level enters into the upper Hubbard band with only slight modifications of the spectral weight. 
Most important, the arc around k = (it, 0) in the upper Hubbard band of the half-filled system is virtually unchanged 
and now forms an electron pocket around k = (tt, 0), as can be seen in panel A. However, for this underdoped system, 
the Fermi level has not yet reached the bottom of the parabolic band around k = (tt/2, it/2). The effect of doping is 
not limited to a rigid shift of the band structure with respect to the half-filled situation. As a matter of fact, some 
new spectral weight is created between lo = —It and lo = — 0.5t at k = (0,0) and between ui — —0.5t and lo = Qt at 
k = (7r/2,7r/2). At higher doping, for x = 0.091 (panel B), this new spectral weight becomes more pronounced and 
starts producing the branches at k = (0,0) — > (it, 0) and k = (0,0) — > (tt/2,it/2) below the Fermi level, which are 
clearly observed in the overdoped system in panel C. 

The results described above provide a deeper insight into the doping process of the t—t' — U Hubbard model, provided 
one assumes that the parameters of the model do not change as a function of doping. Although this assumption is 
widely believed to be appropriate for the doping range considered here, there have been suggestions that the on-site 
repulsion U is constant over a broader doping range. The possibility of a doping-dependent on-site repulsion was 
considered recently on the basis of a spin-density wave (SDW) mean field calculation of the t — t' — t" — U Hubbard 
model to describe the experimental data of the electron-doped system NCCO 0, ■ In contrast to the usual SDW 
calculation, where one self-consistently determines the single particle gap A m f under the assumption that U is a fixed 
parameter, the authors of Ref. |3| consider U as a doping-dependent parameter fixed by the condition that at each 
doping value A m f corresponds to the experimentally-determined pseudogap. In their calculation, which we herewith 
refer to as KMLB, the value of U c s drops sharply upon doping from U e g = Qt at half filling to U c s « it at x ~ 0.15. 
The success of this idea is supported by the fact that the spectral function obtained from this procedure shows 
excellent agreement with the doping evolution of the experimentally observed Fermi surface. The whole scenario, 
however, is based upon the assumption that a varying U e g is indispensable for the reproduction of the experimental 
results. 

For this reason, it is important to verify, whether a doping-dependent U is necessary in order to correctly reproduce 
the spectral features as a function of doping, in particular the doping dependence of the Fermi surface. This is 
conveniently represented experimentally by plotting an intensity plot of the spectral weight at the Fermi energy 
lo = 0. The upper row of Figure |3| (panels A through C) shows the ARPES data taken from 0, whereas the lower 
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FIG. 4: Direct comparison between CPT and KMLB results (see text). The KMLB dispersion is indicated by the solid line. 
In both cases, energies are given in units of t = 1. A: half-filling; B: x = 0.10 (CPT x = 0.091 



row (panels D through F) plots the results obtained by in the CPT calculation. [32} 

Although we only used a minimal set of standard parameters and did not change the parameterization (in particular 
U) as a function of doping, the CPT Fermi surface of the fully correlated Hubbard model qualitatively reproduces 
the experimental result. In particular, we observe electron pockets around k = (jr, 0) in the underdoped region, a FS 
patch at k = (tt/2,w/2) for about x = 0.10 and a large FS centered around k = (w,tt) in the overdoped case. All 
these features are consistent with the experimental results found in |2j as well as in the Hartree-Fock calculation of 
Ref. [3j. As a matter of fact, figure [21 shows that the closing of the gap at higher doping, which is obtained "by hand" 
upon decreasing U within the KMLB, comes about naturally within our numerical treatment of the Hubbard model, 
in which electron correlations are treated more accurately. 

Figure 0| directly compares the CPT and KMLB quasiparticle dispersion at half filling and x = 0.10 (CPT: x — 
0.091). Notice that for the KMLB results only the dispersion is indicated (solid line). In order to have comparable 
energy scales, we reported the energies from |3j in unit of t = 1. 

At half filling (Panel A) both methods show almost identical dispersions for the lowest-energy excitations, since 
both methods have been fitted to the experimental results. However, as a SDW-type mean field method, the KMLB 
results cannot describe the full lower Hubbard band and is exclusively fitted to the low-energy excitations. Hence, 
the lower part of the spectrum between to w — 3i to u> « — At cannot be reproduced within this technique. 

Both techniques describe the closing of the gap as a function of doping. Although the CPT dispersion is much 
weaker and shows a much smaller gap at k = (tt/2, tt/2), the qualitative development as a function of doping is very 
similar. As discussed above, there is additional information provided by CPT in contrast to KMLB, concerning the 
high-binding energy part of the spectrum. In particular, a remainder of the parabolic dispersion at k = (n/2, 7r/2) 
is found at to ~ —5t. This feature does not appear in the mean field calculation, but can be clearly identified in the 
experimental data Q{. 

CONCLUSION 

In conclusion, we have shown that the evolution of the Fermi surface of the electron-doped cuprates is well described 
within the framework of the one-band Hubbard model with nearest- and next-nearest- neighbor hoppings. In particular, 
we have provided indications that a doping-dependent Hubbard repulsion U is not necessary in order to describe the 
doping dependence of the ARPES spectrum. 
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